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In this work, the charm quark mass is obtained from a QCD sum rule analysis of the charmonium system. 
In our investigation we include results from non-relativistic QCD at next-to-next-to-leading order. Using the 
pole mass scheme, we obtain a value of M c = 1.70 ± 0.13 GeV for the charm pole mass. The introduction of a 
potential-subtracted mass leads to an improved scale dependence. The running MS-mass is then determined to 
be m c (m c ) = 1.23 ± 0.09 GeV. 



1. Introduction 

An essential task within modern particle 
physics lies in the determination of quark masses 
which are key parameters of the standard model. 
In the past, QCD moment sum rule analyses have 
been successfully applied for extracting the charm 
and bottom quark masses from experimental data 
on the charmonium and bottonium systems re- 
spectively The fundamental quantity in 
these investigations is the vacuum polarisation 
function Ii(q 2 ): 



(1) 



IW? 2 ) = i J d 4 x e^ x (T{ Jfl (x)M0)}) 

= {%q v - 6v<7 2 )n(<? 2 ), 



where the current is represented by j^{x) = 
(c7 At c)(x). Via the optical theorem, the experi- 
mental cross section cr(e + e~ — > cc) is related to 
the imaginary part of n(s): 



1 cr(e+e 
R ( s ) = 7X2 
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Ql a(e + e — > fi+fi ) 



= 127rimn(s -Me). 



Using a dispersion relation, we can express the 
moments by an integral over the velocity v = 
yl — 4m 2 /s: 
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4m 2 - n(s) 
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The moments can either be calculated theoreti- 
cally, including perturbation theory, Coulomb re- 
summation and non-perturbative contributions, 
or can be obtained from experiment. In this 
way, we can relate the charm quark mass to the 
hadronic properties of the charmonium system. 

A natural choice for the mass appearing in eq. 
(|^) is the pole mass M. In the first part of our nu- 
merical analysis, we will use the pole mass scheme 
to extract the charm pole mass. However, as 
the pole mass suffers from renormalon ambigu- 
ities Q, it can only be determined up to correc- 
tions of order Aqcd- In the second part of our 
analysis we shall therefore use the recently intro- 
duced potential-subtracted mass mps P. From 
this mass definition we can obtain the MS-mass 
more accurately than from the pole mass scheme. 

In the next section, we shall present the pertur- 
bative expansion of the correlator which is known 
up to next-to- next-to-leading order (NNLO). As 
was found in recent analyses of the Upsilon- 
system — |l0|| , the dominant theoretical contribu- 
tions arise from the threshold expansion in the 
framework of non-relativistic QCD (NRQCD). 
These contributions will be derived in section 
3. Afterwards, we shall shortly discuss the non- 
perturbative contributions and the phenomeno- 
logical spectral function. In the numerical analy- 
sis, we shall first explain the method of analysis. 



We will then obtain the pole mass and the MS- 
mass from an analysis in the pole mass and the 
PS-mass scheme respectively. The origin of dif- 
ferent contributions to the error will be carefully 
investigated. We shall conclude with a summary 
and an outlook. 

2. Perturbative expansion 

The perturbative spectral function R(v) can be 
expanded in powers of the strong coupling con- 
stant a s , 

R pt (s) = (s) + ^ J2« (s) + ^ R {2) (*) + ••■ 

From this expression the corresponding moments 
M. n can be calculated via the integral of eq. (g). 
The first two terms are known analytically and 
can for example be found in ref. jy]]. H( 2 \s) 
is still not fully known analytically. However, the 
method of Pade-approximants has been exploited 
to calculate II^ 2 ) numerically, using available re- 
sults at high energies, analytical results for the 
first eight moments and the known threshold be- 
haviour [ p2||lg| ]. This information is sufficient to 
obtain a numerical approximation of I±( 2 )(s) in 
the full energy range. 

The numerical stability of the results can be 
checked in different ways. By choosing different 
Pade-approximants or by selecting a smaller set 
of input data the results for the moments remain 
almost unchanged. Furthermore, some contribu- 
tions to the spectral density like those from inter- 
nal quark loops are known analytically and are in 
very good agreement with the numerical spectral 
density. The above are strong indicators that the 
numerically obtained spectral density comes very 
close to the exact spectral density. 

3. Coulomb resummation 

The perturbative expansion which contains 
terms up to the order 0(a 2 s ) is a good approx- 
imation for high velocities. However, as one 
approaches the threshold region, terms of the 
order v(a s /v) k become increasingly important. 
These terms can be obtained in the framework 
of NRQCD. The correlator is then expressed in 



terms of a Greens function Q: 

N ( 4fc 2 \ 

n(s) = ^ (C h (a s )G(k) + ^Gc(k)) , 

where k = \J M 2 — s/4 and M represents the pole 
mass. The constant Ch(a s ) is a perturbative coef- 
ficient which is needed for the matching between 
the full and the non-relativistic theory. It nat- 
urally depends on the hard scale. The Greens 
function is analytically known up to NNLO || 
and sums up terms of order a™/v n ~ k for n > 
and k = 1, 2, 3. It is crucial for the analysis that 
the result depends on three scales. While the hard 
scale [ihard is responsible for the hard perturba- 
tive processes, the soft scale /x S o/t governs the 
expansion of the Greens function. Furthermore, 
the factorisation scale /i/ ac separates the contri- 
butions of large and small momenta and plays the 
role of an infrared cutoff. 

The Greens function contains two parts: the 
continuum and the poles above and below thresh- 
old respectively. We are interested in both con- 
tributions separately: first, the individual correc- 
tions can be analysed and their error estimated. 
Second, in our numerical analysis we will recon- 
struct the spectral density above threshold and 
we thus need the corresponding spectral density 
at low velocities. In principle, the expressions for 
the energies and decay widths of the poles have 
been calculated. However, in the actual case of 
the charm quark this expansion does not converge 
well. We will therefore choose a different method 
of evaluation. Since the results for the Greens 
function are known analytically, we can evaluate 
their contribution to the moments numerically by 
performing the derivatives. On the other hand, 
by using a dispersion relation, we can obtain the 
continuum from the imaginary part of the corre- 
lator. From the difference we can obtain the pole 
contributions: 

m»- - l|!(4M^)"n(.)| >=o 

_12.(4MTf <fc5^M 

J4M 2 s 

Since we will not evaluate the poles near thresh- 
old, but rather calculate their contributions to 



the moments in a region where perturbation the- 
ory is expected to be valid, the convergence of 
the pole contributions is improved. Nevertheless, 
the poles will give the largest contribution to the 
theoretical moments and thus the dependence on 
the scales will remain relatively strong. The large 
corrections are partly due to the definition of the 
pole mass. These contributions can be reduced 
by using an intermediate mass definition. In this 
analysis we will use the potential-subtracted (PS) 
mass [U where the potential below a separation 
scale ii sep is subtracted: 



5m(fj, sep ) = - 




q|</i«ep 



mpsifJ-sep) = M-Sm(fi sep ). 

As will be seen in the numerical analysis, this 
mass definition leads to an improved scale depen- 
dence and a more precise determination of the 
MS-mass. 

4. Condensate contributions 

The non-perturbative effects on the vacuum 
correlator are parametrised by the condensates. 
The leading correction is the gluon condensate 
contribution which is known up to next-to-leading 
order ]14} |. Furthermore, the dimension 6 and 8 
contributions have been calculated and will be 
included in our analysis p5| , p!q |. It will turn 
out that the condensate contributions are sup- 
pressed when compared to former charmonium 
sum rule analyses and only have little influence 
on the mass. Besides an increase of the theo- 
retical moments from the Coulomb contributions 
we will restrict the moments to n < 7 where 
the non-perturbative contributions are relatively 
small. Since we obtain a larger pole mass than the 
former analyses, the condensates, starting with a 
power of 1/M 4 , are suppressed further. 

5. Phenomenological spectral function 

Experimentally, the first six ^-resonances have 
been observed. Since the widths of the poles are 
very small compared to the masses, the narrow- 
width approximation provides an excellent de- 



scription of these states. To model the contri- 
butions above the 6th resonance, we use the as- 
sumption of quark-hadron-duality and integrate 
over the perturbative spectral density: 

Mn 9x f T k f°° R pt (s) 

(4M2)« a 2 m Q2 2^ E 2n + l+ ^n+l" 

To estimate the continuum contribution we will 
use a threshold in the range of 3.6 GeV < ^/so < 
4.2 GeV with a central value of ^/sq = 3.8 
GeV. However, the most dominant phenomeno- 
logical contributions come from the first two ip- 
resonances resulting in a small influence of the 
continuum even for low values of n. 

6. Numerical analysis 

6.1. Analysing method 

Besides the contributions from the poles of the 
Greens function and the condensates, the theo- 
retical part of the correlator contains the spectral 
density above threshold. Now we will discuss the 
different parts of the spectral density. For high 
velocities the spectral density is well described 
by the perturbative expansion. However, as one 
approaches smaller values of v, the perturbative 
expansion breaks down. The resummed spectral 
density, on the other hand, gives a good descrip- 
tion for low values of v, but becomes unreliable 
for high velocities. For these reasons, we will in- 
troduce a separation velocity v sep w 0.3. Above 
v sep we will use the perturbative spectral den- 
sity. Below v sep , we take the resummed spectral 
density adding the terms which are included in 
perturbation theory but not in resummation. 

In fig. [l] we have displayed the different contri- 
butions. The dotted line represents the pertur- 
bative expansion in LO and NLO. The thin solid 
line also includes the NNLO. Whereas for high 
velocities the perturbative expansion is well con- 
vergent, the importance of the higher corrections 
increases for smaller v. The dashed line is the re- 
summed spectral density and the thick solid line 
the reconstructed spectral density. For the char- 
monium system, there exists a range of intermedi- 
ate values of v where neither the perturbative ex- 
pansion nor the resummation can be trusted. In- 
deed, it can be clearly seen that the reconstructed 




Figure 1. Thick solid line: reconstructed spec- 
tral density; Thin solid line: perturbative spec- 
tral density; dotted line: perturbation theory at 
NLO; dashed line: resummed spectral density. 



Figure 2. Thick solid line: central pole mass 
thin solid lines: M c for fi S oft = 1-05 and 1.5 GeV: 
dashed lines: M c for fif ac = 1.2 and 1.7 GeV: 
dotted lines: M c for fihard = 1.4 and 2.5 GeV. 



spectral density shows a gap at the separation ve- 
locity. To estimate the error we have varied v sep 
between 0.2 and 0.4. The analysis shows that 
though the introduction of the separation veloc- 
ity stabilises the sum rules, the variation only has 
a minor influence on the mass. 

6.2. Pole mass scheme 

Since the perturbative corrections grow for 
large values of n, we will restrict our analysis to 
moments with n < 7. As our analysing method 
needs moments of n > 3 to reconstruct the spec- 
tral density, we will use 3 < n < 7. As the central 
values for our scales we will choose 



fi so ft = 1.2 GeV, 
fihard = 1.7 GeV. 



= 1.45 GeV, 



(3) 



We have set the hard scale equal to the central 
value for the pole mass. For the soft scale, we 
would have liked to choose a somewhat smaller 
value, but the NNLO corrections get out of con- 
trol for fisoft < 1 GeV. The factorisation scale 
lies between the two other scales. In this scheme 
the theoretical moments are dominated by the 
pole contributions. To estimate the error on the 
mass, we have varied the scales within reasonable 
ranges. The result is depicted in fig. [| The thick 



solid line gives the central value for the charm 
mass of M c = 1.70 GeV with the values of eq. 
(Q). The error is dominated by the variation of 
the scales, we obtain 



1.05 GeV < fi soft < 1-5 GeV 
1.2 GeV < fi fac < 1.7 GeV 
1.4 GeV < fihard < 2.5 GeV 



AM C = 100 MeV 
AM C = 50 MeV 
AM C = 40 MeV. 



A significant uncertainty also comes from h-QCD- 
By choosing Aqc d — 330 ± 30 MeV we get an 



Table 1 

Single contributions to the error of M. 



Source 



AM C 



Variation of fisoft 
Variation of fifac 
Variation of fihard 
Threshold so 
Experimental error 
Variation of v sep 
Variation of Aqcd 



100 MeV 
50 MeV 
40 MeV 
10 MeV 
15 MeV 
10 MeV 
50 MeV 



Total error 



130 MeV 



error of AM r = 50 MeV. The results are sum- 



mps 

1.5 r 



Figure 3. Thick solid line: central PS-mass; thin 
solid lines: tops for fX S oft = 1-05 and 1.5 GeV; 
dashed lines: mps for /i/ ac = 1.2 and 1.7 GeV; 
dotted lines: tops for \ihard = 1-4 and 2.5 GeV. 



marised in table [IJ. Adding the errors in quadra- 
ture we obtain the charm pole mass: 



Mr 



1.70 ±0.13 GeV. 



(4) 



Instead of deriving the MS-mass from the pole 
mass, in the following section we will use the PS- 
mass to obtain a more precise value for the MS- 
mass. 

6.3. Potential-subtracted mass scheme 

First, we have to choose a value for the sep- 
aration scale [i sev . This scale should be taken 
large enough to guarantee a perturbative relation 
to the MS-mass. On the other hand, it should 
be smaller than Mv. Both conditions cannot be 
well fulfilled at the same time. As a compromise 
value we will choose fi sep = 1.0 ±0.2 GeV. In this 
scheme the pole contributions from the Greens 
function turn out to be smaller than in the pole 
mass scheme. The contributions from the con- 
densates get more important here and we shall 
restrict our analysis to n < 6 where these correc- 
tions are under good control. 

Using the same central values for the scales (||) 
we obtain mps(fJ- se p = 1-0) = 1.35 GeV and from 
this value a MS-mass of m c (m c ) — 1.23 GeV. The 
introduction of the intermediate mass definition 



leads to a reduced scale dependence: 
1.05 GeV < ^ soft < 1.5 GeV : Am PS = 60 MeV 
1.2 GeV < f ac < 1.7 GeV : Am PS = 40 MeV 
1.4 GeV < fi hard < 2.5 GeV : Am PS = 30 MeV. 
In table | we have listed the individual contri- 

Table 2 

Single contributions to the error of mps and m c . 



Source 



Atops Ato c 



Variation of /z so / t 60 MeV 50 MeV 

Variation of fi fac 40 MeV 40 MeV 

Variation of fi hard 30 MeV 30 MeV 

Variation of fi sep 30 MeV 30 MeV 

Threshold s 30 MeV 30 MeV 

Experimental error 10 MeV 10 MeV 

Condensates 20 MeV 20 MeV 

Variation of v sep 10 MeV 10 MeV 

Variation of A QC d 10 MeV 20 MeV 



Total error 



90 MeV 90 MeV 



butions to the error of mps and m c . Finally, we 
obtain for the masses 



mpsifJ-sep = 1-0) 

m c (m c ) 



= 1.35 ±0.09 GeV 
= 1.23 ±0.09 GeV. 



(5) 

When we use the pole mass from eq. (Q) and 
calculate the MS-mass using the three-loop re- 
lation between these masses |l7||lS|l , we obtain 
to c (to c ) = 1.20 ± 0.17 GeV. This value is in good 
agreement with eq. (|]). This is not self-evident 
as the dominating pole contributions are reduced 
in the PS-scheme and the relative influence of the 
individual contributions is shifted. 

7. Conclusions 

The obtained value for the pole mass lies some- 
what higher than in former sum rule analyses 
Jll],|2C]]. In jl9) the authors used perturbation 
theory to NLO resulting in a value of M c = 
1.46 ± 0.07 GeV. In the second investigation 
PQ| the analysis has been performed in the MS- 
scheme with perturbation theory to NLO. Us- 
ing the NLO relation to the pole mass the au- 
thor obtains m c (m c ) = 1.26 ± 0.05 GeV and 



M c = 1.42 ± 0.03 GeV. The author has also per- 
formed an analysis using resummation in LO with 
a value of M c = 1.45 ± 0.07 GeV. In our analy- 
sis the increased value of the pole mass is essen- 
tially due to large Coulomb contributions which 
have not been included in former analyses. As 
a consequence, the error becomes larger as well. 
In a recent analysis, the pole mass has been es- 
timated from the charmonium ground state at 
NNLO Here the authors obtained a pole 

mass of M c = 1.88±g;f| GeV. This value resulted 
from large corrections of the Coulomb potential 
to the ground state. During the last years, sev- 
eral lattice analyses obtained the following values 
for the MS-mass: 

m c {m c ) = 1.59 ±0.28 GeV [22], 
m c (m c ) = 1.33 ±0.08 GeV §|, 
m c (m c ) = 1.73 ±0.26 GeV ||, 
m c {m c ) = 1.22 ±0.05 GeV §|. 

Whereas the results from |^3|,^5| are in good 
agreement with this analysis, the investigations 
from J22| and (24) obtain higher masses. For the 
time-being, the results are not conclusive and fu- 
ture lattice calculations for the charm mass might 
be of interest. Since other methods reveal sig- 
nificant uncertainties in the determination of the 
quark masses, the sum rules remain one of the 
most precise methods to extract these fundamen- 
tal quantities. 
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